Theory of defect production in nonlinear quench across a quantum critical point 
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We study defect production in a quantum system subjected to a nonlinear power law quench 
which takes it either through a quantum critical or multicritical point or along a quantum critical 
line. We elaborate on our earlier work [D. Sen, K. Sengupta, S. Mondal, Phys. Rev. Lett. 101, 
016806 (2008)] and present a detailed analysis of the scaling of the defect density n with the quench 
rate r and exponent a for each of the above-mentioned cases. We also compute the correlation 
functions for defects generated in nonlinear quenches through a quantum critical point and discuss 
the dependence of the amplitudes of such correlation functions on the exponent a. We discuss 
several experimental systems where these theoretical predictions can be tested. 
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I. INTRODUCTION 

Quantum phase transitions have been widely studied in 
different systems for several years^. Such transitions oc- 
cur when the ground state of a quantum system changes 
due to the variation of some system parameter such as 
pressure 2 -, doping 3 , or magnetic field^. More recently, 
non-equilibrium physics around such critical points has 
also been studied 5 -! 6 -. In particular, quench dynamics 
through quantum critical points has been a subject of 
intense theoretical study in recent years. Such a dynam- 
ics involves the time evolution of a parameter A = A(t) in 
the Hamiltonian of the system which carries it through 
a quantum critical point, characterized by the correla- 
tion length exponent v and the dynamical critical expo- 
nent z, at A = A c . Since the energy gaps between the 
ground and the first excited states vanish at the quan- 
tum critical point, the dynamics of the system necessar- 
ily becomes non-adiabatic in a finite region around this 
point even for an arbitrarily slow quench. This leads 
to the failure of the system to follow the instantaneous 
ground state. As a result defects are produced 7 -^. Most 
of the initial studies of defect production in quench dy- 
namics for various quantum systems have been restricted 
to the case of a linear quench A(t) = Xot/r, where r _1 
is the quench ratoMVVVV^VVV 8 ' 19 ! 20 ' 21 ! 22 . It is 
well known that for a slow linear quench, the defect den- 
sity n ~ T -dv/(zv+i) w j iere d j s th. e dimension of the 
syste m 23 ' . More recently, nonlinear power law quenches 
characterized by X(t) — Xo\t / T\ a sign(t) , where a denotes 
the power law exponent and sign is the signum function, 
have also been studie d 25 ' 26 . In particular, it has been 
shown in Ref. |25| . that if, during the quench, the crit- 
ical point is reached at time t = [t = to ^ 0], then 
the defect density n for such a quench process scales as 



-doif I (azy+l) 



where g 



is a non-universal constant.] 

On the experimental side, trapped ultracold atoms in 
optical lattices have provided ways to realize many in- 
teracting quantum systems with a variety of low temper- 
ature phases separated by quantum critical point o 27 ' 28 . 



These systems provide an easy access to non-equilibrium 
dynamics of its constituent atoms and hence provide ideal 
experimental test beds for quench related studies. Defect 
production has already been studied experimentally for 
a spin-1 Bose condensate 29 . However, a detailed experi- 
mental study of nonlinear quench dynamics has not been 
undertaken so far. 



In this paper we study defect production due to non- 
linear power law quenches in quantum critical systems. 
Our main results are the following. First, we elaborate 
on the work of Ref. [25| and provide a detailed deriva- 
tion of the scaling laws of the defect densities mentioned 
above. Second, we extend the scaling law for defect pro- 
duction through multicritical points, as studied for a lin- 
ear quench in Ref. 21, to nonlinear quenches. Third, 
motivated by the work in Ref. [20L we derive scaling laws 
for defect densities produced during a nonlinear quench 
when the system is taken along a gapless line during the 
quench. Our results extend those in Ref. [2l| and [201 and 
reproduce them as special cases. Fourth, taking the one- 
dimensional Kitaev model as a specific system, we com- 
pute the correlation functions for defects produced dur- 
ing a nonlinear quench. We also provide a general model 
independent discussion of the behavior of such correla- 
tion functions. Finally, we present a detailed discussion 
of possible experimental systems where these theoretical 
results may be tested. 



The organization of the paper is as follows. In Sec. 
[II] we provide detailed derivations for the scaling laws of 
defect density produced during a nonlinear quench. This 
is followed, in Sec. IIII1 by a computation of the defect 
correlation functions. Next, in Sec. IIV1 we provide nu- 
merical studies to corroborate our analytical results. In 
Sec. El we discuss possible experimental systems where 
the scaling laws derived in Sec. [II] can be tested. Finally 
we conclude in Sec. IVTl 
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II. DEFECT PRODUCTION RATE IN A 
NONLINEAR QUENCH 

The density of defects produced in a quench process 
depends crucially on the nature of the phases that the 
system passes through during the quench. Such pro- 
cesses can therefore be broadly classified into three types. 
First, the system may pass from one gapped phase to 
another through an intermediate gapless critical or mul- 
ticritical point. Second, the system may move along a 
gapless critical line in the parameter space so that at 
each point on that line the gap vanishes at a fixed and 
unique momentu m 19 ' 20 . Third, the quench may take the 
system from a gapped phase to another through a gapless 
hypersurface in parameter space as well as in momentum 
spaced. In what follows, we will study defect produc- 
tion during nonlinear quench of the first two types in 
Sees. Ill A[ III Bl and III CI An analogous study for the 
third case, where the system passes through a hyper- 
surface in momentum space, is beyond the scope of the 
present work. 



A. Quench dynamics from one gapped phase to 
another 

We start with the model Hamiltonian for a d- 
dimensional system 

H(t) = £ V4 H s (t) ^, 
fc 

H^t) = (\(t)+b(k))T 3 +A(k)T + +A*(k)T_, (1) 

where b(k) and A(fc) are model-dependent functions, Tj 
denote usual Pauli matrices, \{t) = Ao|t/r| a sign(t) is the 
quench parameter where a — 1 implies linear quench, and 
tpt = (cjgjCjg) represents the fermionic operators. Such 
a Hamiltonian is known to represent several one- and two- 
dimensional spin models such as the Isingi, the XY— and 
the extended Kitaev mode l 18 ' 31 ' 32 . The instantaneous 
energies of the Hamiltonian given by Eq. (fT]) are given 
by 

E(k) = ± ^/(A(0+6(fc)) 2 + |A(fc)| 2 . (2) 

These energy levels touch each other at t = to and fc = fco , 
so that |A(fc)| ~ \k- fc | and \t \ = T|6(fc )/A | 1/a = 
-rg 1 /", where g = |b(fco)/Ao| is a non-universal model- 
dependent parameter. At this point the energy levels 
cross and we have a quantum critical point with v = z = 
1. Note that the critical point is reached at t — only if 
6 (fco) vanishes. 

Let us first consider the case where 6(fco) = so that 
the system passes through the critical point at t = 0. In 
what follows, we shall assume that |A(fc)| ~ |fc — fco| and 
b(k) ~ |fc— fco| Zl at the critical point, where z\ > 1 so that 



E ~ |fc — fco| and z = 1. In the rest of the analysis, we 
set H = 1, and scale t — > tAo, r — > tAo, A(fc) — > A(fc)/Ao, 
and b(k) -» 6(fc)/A . 

We begin by observing that the ground state of the 
system must be (c 1 j, c^£) = (1, 0) [(0, 1)] at the beginning 
[end] of the quench at t = — oo [oo]. Thus the probability 
of defect formation, i.e, the probability for the system 
to be in the excited state at the end of the quench for a 
given state |fc) must be given by 

pz - ^M')l 2 - ( 3 ) 

The density of these defects is thus given by 

where J BZ denotes integration over the Brillouin zone. 

To obtain pj, we study the time evolution of the system 
which is governed by the Schrodinger equation idipz/dt = 
Hpj]^; this leads to the following equations, 

= (|t/r| Q sign(i) + b{k))c l% + A(fc)c 2j j, 
i*2k = -(|t/rrsign(t)+6(fc))c 2fe - + A*(fc)c lfc -, 

(5) 

where we have kept the time dependence of Cjg,^ (t) 

implicit, and ^^^(fc) = dtc^^y To solve these equa- 
tions, we define 

r \ - c p ip dt (\t /T\ a sign{t')+b(k)) 
C lfc _ L lfc ^ 

C = c 2 g e - i / tdt '(l*'/^rsign(t')+fc(fc))^ (g) 

Then substituting Eq. ([6]) in Eq. ([5]) and eliminating c 2 - 
from the resulting equations, we get 

■6 l% - 2z [|t/T| Q sign(t) + b(k)} d 1% 

+ |A(fc)| 2 c 1% = 0. (7) 

Now we scale t — > tT a ^ a+1 ^ so that Eq. (J7J) becomes 

c' lg - 2i [|trsign(i) + 6(fc)W(" +1 )] c' lfe - 

+ |A(fc)| 2 r 2Q /( Q+1 ) c x% = 0. (8) 

From Eq. (|8|) we immediately note that since and c - 
differ only by a phase factor, must be given by 

p % = hm^ 0o |c' lg (t)| 2 = /[Kfc)r^,|A(fe)| 2 r^],(9) 

where / is a function whose analytical form is not known 
for a ^ 1. Nevertheless, we note that for a slow quench 
(large t), becomes appreciable only when the instan- 
taneous energy gap, as obtained from Eq. {2|), becomes 
small at some point of time during the quench. Conse- 
quently, / must vanish when either of its arguments are 
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large: /(oo, a) = f(a, oo) = for any value of a. Thus 
for a slow quench (large t), the defect density n is given 
by 



BZ 



(2n) d 



f[b(k)r^,\A(k)\ 2 r< 



(10) 



and receives its main contribution from values of / near 
k = ko where both b(k) and A(k) vanish. Thus one ob- 
tains, after extending the range of momentum integration 
to oo, 



d d k 



(2tt 

Now scaling k 



\k-h 



o 



'r»+ 1 



|fc-fc | 2 r^T . (11) 



(it - fc )T- Q/(Q+1) , we find that 
^/(|fc|-r^)/^);|fc| 2 ) 



(12) 



where in arriving at the last line, we have used z\ > 1 
and t — > oo . (If z\ = 1, the integral in the first line 
is independent of r, so the scaling argument still holds). 
Note that for a = 1, Eq. (|12[) reduces to its counterpart 
for a linear quench 2 ^. It turns out that the case Z\ < 1 
deserves a detailed discussion which we defer till Sec. 

Next we generalize our results for a critical point with 
arbitrary values of v and z. To this end, we consider 
a generic time-dependent Hamiltonian H\\t] = Hi[X(t)], 
whose states are labeled by \k) and |0) denotes the ground 
state. If there is a second order phase transition, the 
basis states change continuously with time during this 
evolution and can be written as 



(13) 



The defect density can then be obtained in terms of these 
coefficients a%(t) as 



(14) 



Following the analysis in Ref. [23|, one can then obtain an 
expression for the defect density n as 



d d k 



(27T) 



dX(k\ — \0)e 



(15) 



where 5Ej*(X) — E?(X) — Eq(X) are the instantaneous ex- 
citation energies, and we have replaced the sum over k by 
a (/-dimensional momentum integral. We note, following 
Ref. [23|, that near a critical point, 



5E,(X) = AF(A/\k-k \ 



(16) 



where A is the energy gap, z is the dynamical critical 
exponent, and F(x) ~ 1/x for large x. Also, since the 
quench term vanishes at the critical point, A ~ |A| a2y 
for a nonlinear quench, one can write 



5EAX) = \Xr»F'(\X\ azv /\k-k \ z ), 



(17) 



where F'(x) ~ 1/x for large x. Further, one has 
(fc|^|0) = \k-k \~ z G(A/\k-k \ z ) near a critical point, 
where G(0) is a constant. This allows us to write 



A' 



-G'(X az »/\k-k n, (18) 



l> - 1,,-,, 



where G'(0) is a constanti^. Substituting Eqs. (TIT)) and 
(fl8)l in Eq. (fT5|) and changing the integration variables 
to r] = T ai, /( azi,+ ^ \k - ko\ and £ = |ife - fc r 1/(Qiy) A, we 
find that 



(19) 



where C is a non-universal number independent of r. 

Next we focus on the case where the quench term does 
not vanish at the quantum critical point for fc = fco. We 
again consider the Hamiltonian HAt) in Eq. |T]), but now 
assume that the critical point is reached at t = to 7^ 0. 
This renders our previous scaling argument invalid since 

A(fc ) = but b(k ) ^ 0. In this situation, \t /r\ = g 1/a 
so that the energy gap 8E may vanish at the critical 
point for k = ko. We now note that the most important 
contribution to the defect production comes from times 
near to and from momenta near ko- Hence we expand 
the diagonal terms in H^(t) about t = to and k = ko to 
obtain 



H'(t) = 5> f (£) 



afl («-D/»(LJ0)+^)L 7a 



+ A(fc)r+ + A*(fc)r_ ^j(k), (20) 



where b' '(5k) represents all the terms in the expansion of 
b(k) about k = ko, and we have neglected all terms 

R n = (a- n + l)(a-n + 2)..,(a) 

x g^/ a \(t~to)/r\ n s\gn(t)/n\ (21) 

for n > 1 in the expansion of X(t) about to. We shall 
justify neglecting these higher order terms shortly. 

Eq. (|2"0|) describes a linear quench of the system with 
Teff(ct) = t/(q!<?( q;_1 - ) / q ). Hence one can use the well- 
known results of Landau-Zener dynamics 30 to write an 
expression for the defect density, 



d d k 



BZ 



(27T) 



d d k 

BZ W) 



d exp[-7r|A(fc)| 2 r cff (a)](22) 



For a slow quench, the contribution to n comes from k 
near fcg; hence 



M- d ' 2 = (ag^-^/r) 



d/2 



(23) 
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Note that for the special case a = 1, we get back the 



familiar result : 



and the dependence of n on the 



non-universal constant g vanishes. Also, since the quench 
is effectively linear, we can use the results of Ref. H! to 
find the scaling of the defect density when the critical 
point at t = to is characterized by arbitrary v and z, 



(a-l)/a 



/r 



vd/{zv+l) 



(24) 



Next we justify neglecting the higher order terms R n . 
We note that significant contributions to n come at 
times t when the instantaneous energy levels of H'(t) 
in Eq. (|20[) for a given k are close to each other, i.e., 
(i — t )/r ~ A(fc). Also, for a slow quench, the contri- 
bution to the defect density is substantial only when 

is significant, namely, when |A(fc)| 2 ~ l/r c ff(a). Using 
these arguments, we see that 

Rn/Rn-i = (a-n + l)g-^ a (t-t Q )/(nT) 

~ (a-n+l)/(nVr). (25) 

Thus we find that all higher order terms R n >i, which 
were neglected in arriving at Eq. (|23[) . are unimportant 
in the limit of slow quench (large r). 

The scaling relations for the defect density n given by 
Eqs. (|T9|) and (|24|) represent the central results of this 
section. For such power law quenches, unlike their linear 
counterpart, n depends crucially on whether or not the 
quench term vanishes at the critical point. For quenches 
which do not vanish at the critical point, n scales with 
the same exponent as that of a linear quench, but is 
characterized by a modified non-universal effective rate 
7" c ff (a)- If 5 however, the quench term vanishes at the criti- 
cal point, we find that n scales with a novel a-dependent 
exponent advjiazv + 1). For a = 1, T e g(a) = r and 
adv/{azv + 1) = dv j (zv + 1); hence both Eqs. (Til?)) and 
reproduce the well-known defect production law for 
linear quenches as a special case 2 ^. We note that the scal- 
ing of n will show a cross-over between the expressions 
given in Eqs. (fl"§]) and (pM)) near some value of r = r 
which can be found by equating these two expressions; 
this yields To ~ |6(fco)|~ zly ~ 1//Q - For a > 1, the scaling law 
will thus be given by Eq. (O (Eq. (|2Hj)) for r < (>)r . 
We also note here that the results of this section assumes 
that the system passes from one gapped phase to another 
through a critical point and do not apply to quenches 
which take a system along a critical lin o 18 ' 19 . We shall 
deal with this case in Sec. Ill CI 



B. Quench dynamics through a multicritical point 

In this section, we will consider the effect of a nonlinear 
quench in a system of the form given in Eq. {!]), except 
that we now take 

b(k) ~ |jfe- k \ Zl , and A(jfe) ~ \k-k \ Z2 , (26) 



so that the system passes through the critical point at 
t = 0. This will be a generalization of the discussion in 
the first part of Sec. Ill Al where we had Z\ > Z2 with 
Z2 — 1. We will sec below that a separate analysis is 
required if z 2 > Z\. As discussed recently in Ref. [2l|, 
such a condition arises at the multicritical point of a one- 
dimensional spin-1/2 XY model in a transverse field; in 
that model, we find that z\ = 2 and z-i = 3. 

We begin our analysis by comparing the diagonal and 
off-diagonal terms in Eq. (JTJ . From general considera- 
tions, it is clear that defects are mainly produced when 
both |t/r| a sign(t) + b(k) and |A(fe)| are of order 1 or less 
since this is when the instantaneous energy levels given 
by Eq. [2] are close to each other. We now consider the 
forms of b(k) and A(fc) given in Eq. (f2"6"|) . Two possibili- 
ties arise in the limit r — » oo and | A; — /co | ^ . 
(i) If z\ > Z2, then |A(fc)| being of order 1 or less implies 
that b(k) |A(fc)|, namely, b(k) <C 1. In this case, we 
can ignore the term b(k) in Eq. ([T]). This is equivalent to 
saying that the first argument of the scaling function / 
in Eq. © can be set equal to zero. Following arguments 
similar to those leading up to Eq. (fT5|) . we then see that 
the defect density scales as 



-da/[z 2 (a+l)] 



(27) 



which is independent of the value of z\. 
(ii) If Z2 > zi, then |A(fe)| being of order 1 or less im- 
plies that b(k) > |A(fc)|, namely, b(k) > 1. Thus b(k) 
always remains finite as we approach the critical point 
and cannot in general be neglected. In order to have 
|t/r| Q sign(i) + b(k) of order 1 or less, we must there- 
fore have t ^> 1. Let us define a time to as 1 1 / t | Q = 
— sign(to)6(fc) = — sign(to)?7|£; — ko\ Zl , where 77 is an arbi- 
trary non-universal constant. Thus 



|f | = \v\ 1/a r\k-ko\ Zl/a - 



(28) 



In a spirit similar to Eq. (|20[) . we now linearize the func- 
tion |t/r| Q sign(t) + b(k) near t = t , as (|i/r| Q sign(t) — 
|io/r| a sign(t )) = (t- i )a|t /r| Q - 1 /r which, using Eq. 
(|2"B"]). is equal to (t — t )/r c ff (\k\; a) where 

r eff (|fe|;a)=r eff = r^-^lk - k^^-VI" /a. 

(29) 

The effective linearized Hamiltonian can be written as 

H cS = (r 3 (i-io)/rcff + |fc-fc | Z2 Ti) (30) 

and describes a linear quench with r replaced by 
T c ff(|fc|;a). The corresponding defect density is there- 
fore given by the Landau-Zener expression in Eq. (|22p . 
We find that 

P£ - exp [- ttt \k-ko\ [2z2 - Zl(1 - 1/a)] M (Q - 1)/Q /a], 

(31) 
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and 



-da/[2z 2 a+z 1 {l-a)} 



(32) 



Note that the defect density obtained in Eq. (f3"2"| scales 
with an exponent which is independent of the non- 
universal coefficient r\. 

To generalize these results for models with arbitrary 
z\ < Z2 and v, we notice that such models can be de- 
scribed by an effective Hamiltonian H e s(X(t)), where 
A(i) = (t- t )/T cS (\k\;a) and T cS (\k\;a) is given by Eq. 
(f2"9"| . This effective Hamiltonian therefore describes a lin- 
ear quench with a different r e fj for each fc mode and with 
effective dynamical critical exponent Z2 and correlation 
length exponent v. Thus using the arguments of Ref. l23l . 
we get 



d d k 



(27T) 



dA (fc| — |0)e iT "« tfMFa'sW) \ 
dX 

(33) 



where &%(A) ~ |A| Z2,y F'(|A| Z2 7|fc| Z2 ), and F'(ar) ~ 
l/x for large x. Further, one has (fc|^|Q) = 
|fc|~ Z2 G(A/|fc| Z2 ) near a critical point, where G(0) is 
a constant. Using these relations, one obtains n ~ 
/ d d fc / dA'A'^-^AOexp fe|jfc|["(^+i)+*iKi-<*)]/««' 

J A ' dX"X"F(X")j 2 , where A' = A/lfcl 1 /^ and we have 
set 7/ = f without any loss of generality. Then scaling 
\k\ -> |fc| r «"/[a(^"+i)+^i"(i-«)] ) one finally gets 



_—dau /[a{z 2 v+l)+z 1 v{l-a)} 



(34) 



which reduces to Eq. (|3"2"|) for = f. Note that for 
Eq. (|27|) . a generalization to models with arbitrary z^v is 
straightforward, and is given by Eq. (|f 9|) with z replaced 
by ^2- 

Eqs. (|271). dS3) and fjg]) are the main results of this 
section. These results generalize those in Sec. Ill Al to 
defect production for quenches through arbitrary multi- 
critical points. Note that for z\ — Z2, Eq. (|2"Tj) and Eq. 
([32)) agree for any value of a, giving n ~ T - dQ /[ 2 2("+i)] 
Further, for the case of linear quenching, a — 1, these 
equations agree for any value of Z\ and 22, giving n ~ 
T -d/(2z 2 ) wn j cn nas been recently obtained in Ref. l2~il . 



C. Quench dynamics along a gapless line 

Recently quench dynamics in a one-dimensional XY 
model in the presence of a spatially modulated transverse 
magnetic field has been studied in Ref. [2(1 Such a model 
is described by the Hamiltonian 



^+0 



where J and 7 are respectively the strength of and the 
anisotropy in the nearest neighbor spin-spin interactions, 
a a (a = x, y, z) denote the Pauli matrices, and h and S 
denote the uniform and alternating components of the 
magnetic fields respectively. The phase diagram of this 
model is discussed in detail in Ref. [20- It was pointed out 
that quenching the anisotropy parameter j(t) = 7oi/r 
linearly while sitting at the paramagnetic phase deter- 
mined by the condition h 2 = 8 2 + J 2 , leads to a time 
evolution of the system along a gapless line. It was also 
shown, via mapping this spin model to a system of Majo- 
rana fermions by a Jordan- Wigner transformation, that 
the evolution of the model described by Eq. (|35j) can be 
represented by an effective Hamiltonian given by^ 

H eG (k;t) = X>l(W"3 + Jfc 2 n)Vfc, (36) 

k 

where ipk = {cik,c-2k) is the usual two component 
fermionic field, j(t) = ^(t)J/VS 2 + J 2 and J = 
J 2 /y/d 2 + J 2 . The quench dynamics of this model was 
studied in Ref. for the linear quench j(t) = jot/r us- 
ing the Landau-Zener formalism. It was found that the 
defect density scales as 



--1/3 



(37) 



Note that since for this model z = v = 1, the naive 
expectation according to the analysis of Sec. Ill Al is to 
have n ~ 1/y/r. This result therefore clearly points out 
the necessity of extending the analysis of Sec. Ill Al for 
quenches along gapless lines in parameter space. In what 
follows, we shall only restrict ourselves to quenches where 
the gap vanishes at the same momentum value fco ; fco = 
for the present case. 

We start with a generic d-dimcnsional model described 
by a Hamiltonian 

H cS (k;t) = Y / 4( X (t)\k\ a T 3 + A Q \k\ b T 1 ^ k ,(38) 

k 

where X(t) = Ao|i/V| a sign(t) is the quench parameter, a 
and b are arbitrary exponents, and we have taken fco = 
for clarity. Note that d = a = 1 and 6=2 corre- 
sponds to the XY model studied in Ref. [2(| while a = 
and d = b = 1 corresponds to the one-dimensional Ki- 
taev model studied in Ref. 18. For a / 0, the system 
passes along a gapless line during the quench. We study 
the time evolution of the model in a manner similar to 
that described in Sec. Ill Al After some straightforward 
algebra, one obtains the equation for the evolution of 
c ifcW = c ifeW ex P(* J t ^>7^rsign(t')|fc| a ) as 

■6 l% - 2i\t/r\ a sign(t)\k\ a + A 2 \k\ 2b = 0. (39) 

Next we define tv = r/fc a / a and scale t — > tr^ +1 ' in 
Eq. ([39j) to obtain 



(a*a* +1 -ay*« +1 )-(h-(-iy5)a* (35) 



life 



2i|t| Q sign(i) cf 1& 
\&\l \k\ 



(40) 
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From Eq. (|40)). we find that the probability of defect 
formation for a given momentum k must be given by 

Pi = hm |c',(t)| 2 = /[A^i^a/^!)^/^ 

(41) 

where /[oo] = 0. The defect density therefore becomes 
d d k 



DZ 



^/(A^lfcr^r). (42) 



Using the same logic as outlined in Sec. Ill Al we scale 
|fc| _> T b {a +i)-a |fc| an( j g e t 



n ^ t i>c<»+i)-° 



(43) 



This result generates the scaling of the defect density 
derived in Ref. (n ~ r - d /( 2b - a )) for the special case 
a = 1, and that of the one-dimensional Kitaev model 
studied in Ref. [Hj for a = b = z = 1 and a = 0. 

Finally, we generalize the result in Eq. (|4"3"|) to sys- 
tems where the energy difference between the ground 
and excited states vanishes along the gapless line as 
A = A ac |fc| a . Note that for the quenches treated in Sec. 
Ill Al c = zv. Here, however, since the quench takes place 
along a gapless line, c need not have the same interpreta- 
tion and can be system specific. Exactly at the quantum 
critical point A = A c , the energy gap vanishes as A ~ |fc| b . 
Thus b is to be interpreted as the dynamical scaling expo- 
nent z in the present case. Then using the same scaling 
argument as in Sec. Ill Al we can express the defect den- 
sity n using Eq. (|15p . However, in the present case the 
energy gap scales as 



E k (X)-E (X) = AF(A/\k\ b ), 

(fc|^|0) = l/|fc| b G(A/|fc| b ), 



(44) 



where F(x) — 1/x for large x, and G(0) is a constant. 
Using Eq. (|33)). we see that 



E k {X)-E (X) 



= A ac |fc|°F(A ac |fc|7|fc| b ), 

A ac-1 / A ac|||a^ 



|fc|(6-a) 



(45) 



Substituting these in Eq. (fl"5|) . the defect density pro- 
duced in this system is found to be 



d d k 
(2^p 



A" 



-G 



X ac 



k\{b-a) y\k\b-o 
ir \- x d\'\' ac \%\ a F(\"> c \%\° /\Y'-> , 4i)j 



Defining new variables £ = A|fc|( a b )/ ac and \k \ 
|fc|r ac/[(& - Q)+bQcl , we get 



n ~ r 



- dac / [(b— a)+bctc] 



Eq. (|47|) is one of the central results of this work, and 
it generates all the previous scaling laws for both linear 
and nonlinear quenches through critical lines and points 
(but not through multicritical points) as special cases. 
For a = 0, b = z and c = zv, we recover the scaling law 
Eq. ()19j) for a nonlinear quench, whereas for a = c = 1, 
we obtain the scaling law derived for a linear quench in 
Ref. US- 



UI. DEFECT CORRELATION FUNCTIONS 

For the purpose of computation of defect correlation 
functions, we are going to restrict ourselves to the class of 
d-dimensional models given by H(t) in Eq. ([T]). As men- 
tioned before, many standard spin models in one and two 
dimensions can be mapped, via standard Jordan- Wigner 
transformations^, to such fermionic models described by 
H(t). Let us denote the ground and the excited states 
of H{t) before the quench (at t — — oo) by |0)g and 

respectively for a given value of k. Then the state of the 
system after the quench (at t = oo) is given byi£ 



(48) 



Using this, one can compute the defect correlation func- 
tions for these models. These correlation functions are of 
two types. They can either vanish at the origin, as in the 
case of the two-dimensional extended Kitaev models, or 
can be written a a 10 ' 18 

(CV) = -5? i0 +C [ d d kf[\A(k)\ 2 T^]g(k-f), (49) 
Jo 

where Op — iipftipfi+p, ipa denotes the field operators for 
Majorana fermions, g(k • r) is a system specific function 
independent of t, C denotes a system and dimension spe- 
cific constant which will be unimportant for subsequent 
discussions, and we have used Eq. §§§ to obtain the value 
of Since for a slow quench, is appreciable only near 

k = ko, we expand A(fc) about ko, scale the momentum 
components k[ = (k — ko)iT a ^ a+1 \ and extend the range 
of integration to oo to get 



(Of) = —Sffi + 



c 



-da/(ct+l) 



d d k'f[\k\']g(k -f), 



(50) 



(47) 



where = ri/T a ^ a+1 \ Thus we find that quite gen- 
erally, for the class of models whose defect correlation 
functions do not vanish at the origin, 

drv 

ln(l + (O ?=0 » = ln(C')-— -Mr), (51) 

a + 1 

i.e., the logarithm of the deviation of the amplitude of 
these correlation functions at the origin from —1 is a 
linear function of ln(r) with a slope of —da/ (a + 1). 
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r 

FIG. 1: Plot of (O r ) vs r for a = 2 (black dot-dashed line), 
a = 3 (red dashed line), a = 4 (blue solid line), and r = 20. 

We now compute the correlation function for a specific 
model, namely, the one-dimensional Kitaev model [I8II25I 
[33l | which has the Hamiltonian 

H = ( J i S i S f+i + J ^ S t S ti) , (52) 

iGeven 

where J\ and J2 denote the nearest neighbor interaction 
strengths, and Si denotes the spin at site i. Using the 
standard Jordan- Wigner transformation, this Hamilto- 
nian can be mapped on to a free fermionic Hamiltonian 
0111 

H = ^2 ip' k H k ip k , where 

k 

H k = -2 (J_ sin(fc) r 3 + J+ cos(fc) r 2 ). (53) 

Here J± = J\ ± J2, and ipk = (ci(fc), 02(h)) are the 
fermionic fields. The Hamiltonian is changed in time by 
varying the parameter J_ keeping J + fixed. The defect 
correlation function for this model is given by [T|| 

2 f* 

(O r ) = — S r o + — / dh Pk cos(fcr). (54) 

7T Jo 

Thus we find that the defect correlation functions have 
the same form as in Eq. (j49| with C — 2/tt and g = 
cos(fcr). A plot of the correlation function as a function 
of r, sans the <5-function peak at the origin, is shown in 
Fig. [T]for t = 20 and several representative values of a. 



IV. NUMERICAL EVALUATION OF DEFECT 
DENSITIES 

In this section, we provide numerical studies of the one- 
dimensional Ising and Kitaev models to supplement our 




1 ■ ' 

0.2 X 1 1 ■ 1 1 r 

14 8 12 

a 

FIG. 2: Variation of the defect density n with the quench 
exponent a for representative values of r = 10 (black solid 
line), t — 15 (red dashed line) and r = 20 (blue dotted line). 
A polynomial fit of the form n = aa b yields exponents which 
are very close to the theoretical result 1/2 for all values of r 
(see text for details). 

analytical results. First we consider the one-dimensional 
Ising model in a transverse field described by 

Rising = - J C£ S*S? +1 + g Sf), (55) 

i i 

where J is the strength of the nearest neighbor interac- 
tion, and g — h/ J is the dimensionless transverse field. 
In what follows, we shall quench the transverse field as 
g(t) = |i/r| Q sign(t) and compute the density of the re- 
sultant defects. 

We begin by mapping Hi slng to a system of free 
fermions via a standard Jordan- Wigner transformation 1 

H' = - ^((ff-C08(fe))7S + Bin(fe)Ti)V> fc . (56) 

k 

If the external magnetic field g is varied with time as 
g(t) = go\t/T\ a sign(t), then the system will go through 
two quantum critical points at g — 1 and — 1. The energy 
gap vanishes at these quantum critical points at h = 
ho = and ir. As a result, defects are produced in non- 
adiabatic regions near these points. For this model, the 
quantum critical point is at t = t ^ and z = v = 1. 
Hence, T G fr — r/a for both the quantum critical points. 
From Eq. |2~4"1) . therefore, we expect the defect density 
produced in this system to be given by n ~ (r/a)^ 1 ^ 2 . 

To verify this expectation, we numerically solve the 
Schrodinger equation idtifik = Hk(t)ipk and obtain the 
probability pk for the system to be in the excited state. 
Finally, integrating over all k within the Brillouin zone, 
we obtain the defect density n for different values of a > 
1 with fixed r. The plot of n as a function of a for 
t = 10, 15 and 20 is shown in Fig. [5J A fit to these 
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-0.667 
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FIG. 3: Plot of ln(n) vs ln(r) for the one- dimensional Kitaev 
model for a — 2 (black solid line), a — 4 (red dotted line), 
a = 6 (blue dashed line) and a = 8 (green dash-dotted line). 
The slopes of these lines agree reasonably with the predicted 
theoretical values —a/ (a + 1) as shown in the table. 



curves gives the values of the exponents to be 0.506 ± 
0.006, 0.504 ± 0.004 and 0.505 ± 0.002 for r = 10, 15 
and 20 respectively which are remarkably close to the 
theoretical value 1/2. The systematic positive deviation 
of the exponents from the theoretical value 1/2 comes 
from the contribution of the higher order terms neglected 
in the derivation of Eqs. (j2"3")l and (j2~4"|) . We note that the 
region of validity of our linear expansion, as can be seen 
from Fig. [21 grows with r which is in accordance with 
the result in Eq. (j25j) . 

Next, we consider the one-dimensional Kitaev model 
which is governed by the Hamiltonian in Eq. (|52p . As 
mentioned in Sec. 1541 such a model can also be mapped 
on to the free fermionic Hamiltonian given by Eq. (|53[) . 
This system passes through the quantum critical point 
at J_ = for k = n/2 when J_(t) = J_ |i/r| Q sign(£) is 
varied nonlinearly with time. Here the quantum critical 
point is at t = 0. Thus from Eq. (fTTj)) we expect the 
defect density n ~ T -«/( Q + 1 ) since v = z = 1 for this 
system. To check this prediction, we numerically solve 
the Schrodinger equation id t ip{k) = H'(k;t)ip(k,t) and 
compute the defect density n — J Q dk/irpk as a function 
of the quench rate r for different a with fixed J + / J = 1. 
A plot of ln(n) vs ln(r) for different values of a is shown 
in Fig. [3] The slope of these lines, as can be seen from 
Fig. [3l changes from —0.67 towards —1 as a increases 
from 2 towards larger values. This behavior is consistent 
with the prediction of Eq. (fTT)| . The slopes of these lines 
also show excellent agreement with Eq. (|19[) as shown in 
the inset of Fig. [3] 

Finally, we illustrate the expressions in Eqs. ([2T| and 
P2p by taking two one-dimensional models governed by 
Eqs. dXD and with z x = 2, z 2 = 1 and zy = 1, z 2 = 2 
respectively. Setting a = 4, we numerically carry out 
the time evolutions for different values of the momentum 



-0.5 



-1.0 



f.1.5 



-2.0 



-2.5 



3.0 3.1 3.2 3.3 3.4 
ln(t) 



FIG. 4: Plots of ln(n) vs ln(r) for models with z\ — 2 and 
Z2 = 1 (red dotted lower line) and Z\ = 1 and 22 = 2 (black 
solid upper line), for d = 1 and a = 4. The slopes of the 
lower and upper lines are —0.828 and —0.301 which compare 
reasonably with the predicted theoretical values of —4/5 = 
-0.8 and -4/13 = -0.308 respectively. 



k and then integrate to compute the defect density as a 
function of r. The results are shown in Fig. 21 reasonable 
agreement is obtained with the theoretical values of the 
exponents given in Eqs. (j27| and (|32|) . 



V. EXPERIMENTS 

The generality of our results allows for their verifica- 
tion in several realizable experimental systems. We note 
that all our results have been obtained at zero temper- 
ature with the assumption that the system does not re- 
lax significantly during the quench process and till the 
measurement of the defect density has been performed. 
This might seem too restrictive. However, we would like 
to point out that systems of ultracold atoms in opti- 
cal or magnetic traps and/or optical lattices can eas- 
ily satisfy the required criteria since they have a very 
long relaxation time which often gets close to the sys- 
tem lifetime^!. We list some possible experiments briefly 
here. First, there has been a concrete proposal for the 
realization of the Kitaev model using an optical lattice 2 ^. 
In such a realization, all the couplings can be indepen- 
dently tuned using separate microwave radiations. In 
the proposed experiment, one needs to keep J3 = 
and vary Ji( 2 ) — J(l ± |i/r| Q sign(t))/2, so that J + re- 
mains constant while J_ varies in time. The variation 
of the defect density, which in the experimental set-up 
would correspond to the bosons being in the wrong spin 
state, would then show the theoretically predicted power 
law behavior in Eq. (|19[) . Secondly, a similar quench 
experiment can be carried out with spin-1 bosons in 
a magnetic field described by an effective Hamiltonian 
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HoS = c 2 n (S) 2 +ciB 2 (S 2 ) — , where c 2 < and n is the 
boson density. Such a system undergoes a quantum phase 
transition from a ferrom agnetic state to a polar conden- 
sate at B* = \/ |c2 |«.o / c i ■ A quench of the magnetic field 
B 2 = B 2 \t/r\ a would thus lead to a scaling of the defect 
density with an effective rate r e fj(a) = r/ (aj' a_1 '/°), 
where 5 = Ic^l^o/ci. A measurement of the dependence 
of the defect density n on a should therefore serve as a 
test of the prediction in Eq. (f24|) . Finally, spin gap dimer 
compounds such as BaCuSi20g are known to undergo a 
singlet-triplet quantum phase transition at B c ~ 23. 5T 
which is known to be very well described by the mean- 
field exponents z — 2 and v = 2/3 34 . Thus a nonlinear 
quench of the magnetic field through its critical value 
B = B c + B \t/T\ a sign(t) should lead to a scaling of the 
defects n ~ T ~Sa/(4a+3) m ^ — 3 T n the experiment, 
the defect density would correspond to residual singlets 
in the final state which can be computed by measuring 
the total magnetization of the system immediately after 
the quench. We note that for these dimer systems, it will 
be necessary to take special care to achieve the criterion 
of long relaxation time mentioned earlier. 



VI. CONCLUSIONS 

In conclusion, we have studied defect production in 
quantum critical systems for an arbitrary nonlinear 
power law quench. We have shown that the defect pro- 
duction rate depends crucially on whether the system 
passes from one gapped phase to another or along a crit- 
ical gapless line during the quench. We have obtained 
general scaling laws for defect densities produced during 
the quench for both these cases, and have verified these 
laws by numerical studies of one-dimensional systems. 
We have also computed the defect correlation functions 
for a class of <i-dimensional models and have discussed 
the scaling of the amplitude of these functions with the 
quenching rate. Finally, we have discussed several exper- 
imental systems where these results can, in principle, be 
tested. 
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